Interface induced giant magnetoelectric coupling in multiferroic superlattices 
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The electric and magnetic properties of (BaTi03)n/(CaMn03),i short-period superlattices are 
studied by the first-principles calculations. The local electric polarizations in the CaMnOs layers 
are significant, comparable to that in the BaTiOa layers. Remarkably, the electric polarization is 
almost doubled when the spin configuration changes from antiferromagnetic to ferromagnetic in the 
superlattices, indicating a giant magnetoelectric coupling. This enhancement of the magnetoelectric 
coupling is due to the suppression of the antiferrodistortive mode in the CaMnOs layers at the 
interfaces. 

PACS numbers: 77.22.-d, 77.22.Ej, 77.80.-e, 77.84.Lf 
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The emergence of new phenomena at artificial het- 
erostructure is currently at the center of scientific and 
technological interest [l|45|. Because a large variety of 
degrees of freedom such as spin, charge, structural or- 
derings can be found in ABO3 perovskite, it presents an 
ideal playground to explore the interactions among var- 
ious orderings that could potentially lead to enhanced 
functionalities, and multifunctional materials [ll|5|. 

Multiferroics are among these multifunctional ma- 
terials that have attracted intensive interests re- 
cently [6|]. The ferroelectric/ (anti)-ferromagnetic super- 
lattices (SLs) may have strong electric polarization, mag- 
netic ordering and magnetoelectric coupling simultane- 
ously which is rarely seen in single phase bulk materials 
because of the symmetry restrictions f7|. Convention- 
ally, the magnetoelectric coupling can be introduced by 
the mechanical boundary condition through the in-plane 
strain or the electric boundary condition through charge 
continuity satisfied by the constituent materials [7] . How- 
ever, the interfaces separating different constituent ma- 
terials can carry distinct instabilities originated from in- 
dividual parent bulk materials [l| . We show that the in- 
terfacial competition of instabilities [8] may lead to sur- 
prisingly enhancement of the magnetoelectric coupling, 
which arising from the interfacial atomistic effects, can 
be much larger than that through the continuum media 
coupling. 

AMnOs (A=Ca, Sr, ...) are good candidates as build- 
ing blocks for the multiferroic supperlattices, because 
they have several competing instabilities [9| coupled to 
the magnetic ordering. In AMnOa, Mn atom has par- 
tially occupied d-orbitals, and the Mn-O-Mn bond angle 
is crucial in controlling the magnetic interactions, which 
at the same time also strongly couples to the oxygen octa- 
hedral rotation (AFD) and ferroelectric (FE) soft mode. 
The coupling between the spin, the AFD modes and the 
FE modes depends on the relative energetics of these 



instabilities. Unfortunately in bulk AMnOs, there is a 
strong AFD instability associated with a large oxygen 
octahedral rotation that suppresses the FE mode [91]. As 
a result, the Zmear magnetoelectric (spin-ferroelectricity) 
coupling is usually found to be weak. 

In this letter, we demonstrate, through first-principles 
calculations, that the energetics of the AFD, FE and 
magnetic ordering can be drastically modified [1(| by in- 
terface engineering to favorite the magnetoelectric cou- 
pling. We take the CaMn03(CMO)/BaTi03(BTO) SLs 
as our model systems. We find that the MnOg octa- 
hedral rotation will be strongly suppressed by the neigh- 
boring BaO layers, leading to the enormous enhancement 
of magnetoelectric coupling as well as FE at interfaces. 
This enhancement will be strengthened with the increas- 
ing density of interfaces and reaches its maximum at 
the shortest SL, (i.e., n—\) where one observes a huge 
change of electric polarization between the antiferromag- 
netic (AFM) and ferromagnetic (FM) states. To better 
clarify the mechanism for its underlying physics, an ef- 
fective Hamiltonian model is also developed. The model 
explicitly shows that an increased magnetoelectric cou- 
pling develops with a completely suppressed AFD, which 
consistently explains our computational results. 

The calculations are based on the standard density 
functional (DFT) theory with spin-polarized generalized 
gradient approximation, implemented in the Vienna ab 
initio simulations package (VASP) [ll|, [l^. We adopt 
the Perdew-Burke-Ernzerhof functional revised for solids 
(PBEsol)[3]. The on-site Coulomb C/=4.0 eV and ex- 
change interaction J =0.88 eV are used for the Mn 3d 
electrons [IJI . We use the plane- wave basis and projector 
augmented- wave pseudopotentials [15| . A 500 eV energy 
cutoff and 6x6x2 k-point meshes converge the results 
very well. All ionic coordinates are fully relaxed until 
the Hellman-Feynman forces on the ions are less than 1 
meV/A. The electric polarizations are computed using 
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the Berry phase theory of polarization 

The calculated lattice constants of cubic CaMnOs 
and BaTiOs are 3.73lA and 3.987A respectively using 
PBEsol functional |13| , compared with the experimental 
values 3. 73 A and 3. 993 A, which are much more accurate 
than the those obtained from LDA. This is very impor- 
tant for studying ferroelectrics, whose properties are very 
sensitive to the lattice constants. 

The ground state of orthorhombic (bulk) CaMnOa is 
G-type antiferromagnetic. We find a small FE instabil- 
ity related to an unstable polar mode in the high symme- 
try cubic phase, with an imaginary frequency Wi^£;=3.43i 
cm~^ and a much larger antiferro-distortive (AFD) insta- 
bility associated to a nonpolar oxygen rotational mode 
with an imaginary frequency a;AF_D=219i cm^^, consis- 
tent with previous calculations [9|. As a result, the octa- 
hedra rotate a angle about 8.94°. In contrast, BaTi03 is 
highly resistant to oxygen octahedral rotations and ex- 
hibits a robust FE state at room temperature. 

The lattice mismatch between CMO and BTO is 
about 6.59%, which might be too large to grow 
high quality CMO/BTO SLs. One may grow the 
CMO/BTO SLs on the NdGaOs substrate to reduce 
the lattice mismatch between CMO and BTO to about 
3.37%. We therefore fix the in-plane lattice con- 
stants of the SL to those of NdGaOa substrates (3.86 
A) in the calculations. We have studied two short- 
period SLs structures (BaTi03)i/(CaMn03)i (1:1) and 
(BaTi03)2/(CaMn03)2 (2:2). In the calculations, we fix 
the symmetry of the SLs to the space group P4bm, allow- 
ing the MnOe and TiOg octahedra to rotate. We neglect 
the tilting of MnOe octahedra (rotations about an in- 
plane axis), because oxygen tilting requires a coherent 
pattern of tilts that would propagate into the BTO unit 
cells, where octahedral rotations are unfavorable. 

To determine the ground-state spin structure, we cal- 
culate the total energies of a set of spin configurations 
(SCs) for the 1:1 and 2:2 SLs with full relaxations of the 
(electronic and lattice) structures. The results show that 
the stablest SCs arc G-typc AFM for 1:1, and A- type 
AFM for 2:2 SLs. The calculated magnetic moments of 
Mn ions are about 2.9 /is, whereas Ti ions have negligible 
induced magnetic moments. We further fit the exchange 
integrals to a Heisenberg model H = 1/2 J2i , Jij^i ■ Sj, 
assuming nearest neighbor coupling between the Mn ions. 
In the 1:1 SL, the intra-layer exchange interaction Jintra 
— -4.8 meV//i^, and the inter-layer exchange interac- 
tion Jintcr — -0.86 meV/^g. In the 2:2 SL, the ex- 
change interaction in the Mn02 layer that sandwiched 
by the two CaO layer J^^^^ = 8.61 meV/^g, and in the 



Mn02 layer that adjacent to the BaO layer Jj 



(2) 



-1.11 



r(l) 



The interlayer exchange interaction J^^^^j. 



meV//i|. 

-24.54 meV/fj,% and Jj^tcr =-0.04 meV//i| respectively. 
Here, the magnetic moments S are normalized to 1. The 
^ntra ^'^ ^he 2:2 SL is shghtly frustrated. The Neel tem- 
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FIG. 1: (Color oline) (a) The local electric polarizations and 
(b) the AFD rotations of the Mn02 and Ti02 layers in the 
1:1 SL of the G-type AFM (red line) and FM states (blue 
line), (c) The local electric polarizations and (d) the AFD 
rotations of the Mn02 and Ti02 layers in the 2:2 SL of the 
G-type AFM (red line) and A- type AFM states (blue line). 



peratures of the 1:1 and 2:2 SLs are 52 K and 59 K re- 
spectively calculated from Monte Carlo simulations |l7|. 
This is because the Neel temperature of bulk CMO is low 
(^ 130 K) 18]. The Neel temperature of the SLs can be 
enhanced by using other AMn03 compounds with higher 
Neel temperatures. 

Figure[T](b), (d), depict the local AFD associated with 
the TiOg or MnOg octahedral rotation in 1:1 and 2:2 SLs 
respectively. Clearly, the MnOe octahedral that sand- 
wiched between two CaO layers in 2:2 SL keeps a large 
rotation angle around 12° and TiOe octahedra that sand- 
wiched between two BaO layers has a rotation almost 
around 0°. This is in consistent with the fact that CMO 
bulk has a strong AFD instability and BTO strongly re- 
sists the octahedra rotation. What is more intriguing, 
however, is the behavior of the interfacial layers of MnOe 
octahedra. Because of the presence of BaO on one side 
and CaO on the other side, MnOg at interface is exposed 
to a strongly broken mirror symmetry. As a result, this 
MnOg octahedra is found to have significantly reduced 
rotation angle which is around 4.6° and 3.0° in 2:2 and 
1:1 SL respectively. 

In bidk CaMn03, it is the strong AFD that suppresses 
the FE instability from condensation at ground states. 
Because of the greatly reduced interfacial MnOg octa- 
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FIG. 2: (Color oline) The total energies as functions of the 
Mn02 intra-layer exchange energy Jintra for (a) bulk CaMnOs 
and (b) 1:1 SL through AFD rotation and FE modes. We also 
show the AFD rotation angle, and the polarization of the FE 
mode in the 1:1 SL. 



hedra rotation, together with strong polarization in the 
BaTiOa layers, one might expect that the FE can develop 
in the CaMnOa layer of the SL. Indeed, we find that both 
1:1 and 2:2 SLs become multiferroic with spontaneous 
polarizations of 21.30 /iC/cm^ and 27.52 fiC/cm^ respec- 
tively computed by Berry phase formalism. To gain more 
insight, we calculate the layer polarization along [001] di- 
rection pz = ^j {e/n)Z*Xi based on the linear approx- 
imation involving effective charges Z* and small ionic 
distortions Xi of each atoms i in the cell from a higher 
central symmetric nonpolar reference structure. The ef- 
fective charges are obtained from the cubic CaMnOa and 
BaTiOs phase by finite difference method. The resulting 
layer polarizations are presented in Fig.l. It can be seen 
that Mn02 centered layers becomes polarized to almost 
the same degree of Ti02-centered layers in both 1:1 and 
2:2 SLs. 

The development of multiferroic behavior in the 
CaMnOa component in the sperlattices has now been es- 
tablished. It provides us an opportunity to study the 
magnetoelectric coupling directly. To this end, we com- 
pare the electric polarizations under various spin con- 
figurations. As one of the most interesting results, we 
find that the electric polarization undergoes a huge in- 
crease from about 21 ^C/cm^ to about 38 /iC/cm^ in 
1:1 SL when the Mn02 intra-layer spin configuration is 
changed from antiferromatic (C-type or G-type AFM) to 
ferromagnetic (FM or A- AFM). In 2:2 SL, we observe a 
similar however weaker enhancement of polarization from 
about 19 fiC/cin^ to about 30 /iC/cm^. The above result 
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FIG. 3: Schematic shows of the energy surfaces of the effective 
Hamiltonian Eq. [T] in the case of (a) the coupling coefficient 
7 < 7c, and (b) the coupling coefficient 7 > 7c, when both 
a <0, and /3 <0. 



indicates that the magnetoelectric coupling is giant! 

At a closer inspection of what happen locally in the 
SLs, the above drastic polarization change is also found to 
be accompanied by a almost homogeneous increased layer 
polarization for both Mn02-centered and Ti02-centered 
layers in both 1:1 and 2:2 SLs and a further drop of the 
interfacial MnOg octhedra rotations. 

It is well known that the Mn-O-Mn exchange angle 
is essential in controlling the magnetic ordering |19[. In 
multiferroics where spin, FE and AFD instabilities are 
all present, Mn-O-Mn angle can be strongly coupled to 
MnOg octahedral rotation and FE soft phonon modes. In 
Fig. [UJa), (b), we present the change of the total energy 
as a function of the intra-layer exchange integral Jintra in 
bulk CMO and 1:1 SL respectively. It can be seen that 
intra-layer spin ordering can be continuously tuned from 
AFM( Jintra < 0) to FM( Jjntra > 0) by increasing either 
the MnOe octahedral rotation (red) or the amplitude of 
the FE soft phonon mode in both bulk CaMnOs and 1:1 
SL. In bulk CMO, the intra-layer AFM-FM phase transi- 
tion is realized by the increased octahedral rotation which 
has a much lower total energy than that of increased FE 
soft phonon mode. However in 1:1 SLs, the MnOg octa- 
hedral rotation is lagerly suppressed by the interface be- 
cause it becomes energetically unfavorable [see Fig.[21[b)]. 
Thus, AFM-FM phase transition can be only driven by 
the increased FE soft phonon mode which results in the 
observed enhanced magnetoferroelectric coupling. 

It is clear that the enhancement of the magnetoelectric 
coupling originates from the interface where the relative 
energetics of spin, FE and AFD are changed compared 
with bulk CaMnOa. The AFD instability is largely sup- 
pressed thus the spin-polarization coupling strength is in- 



creased. To further understand how the interface tunes 
the competetion between the AFD and FE modes, we 
resort to an effective Hamiltonian. For bulk CaMnOa in 
a uniform phase, the effective Hamiltonian is, 



^({ur,u^}) = E„~Y.J^i(^)^^-^i 



(1) 

+ aul + I3ul + ?7< + KU^ + "fulul , 

where, ur and u,^ are the FE and AFD modes re- 
spectively. The phonon frequencies a — ^mrto'^ — 

Sy W^^^ ■ Sj and /3 = ^m^ujl - J2^J W^^^ ' ^j 
strongly depend on the spin configurations of the sys- 
tem. The anharmonic terms k, rj are positive. 7 
which describes the coupling between the FE and AFD 
modes, is also positive. At the energy minima, we have 
dE{{ur,u^)/dur=0, and dE{{uT,u^)/du^^O, i.e.. 



aur + 277Up + 7Uru = 



/3u, 



V -<- 2ku^ 



7U^Ur = 



(2) 
(3) 



It is easy to see from the equations that if a > 
(/3 >0), we have ur=0 (u^=0). The interesting cases 
are that both a and /3 are negative. When the cou- 
pling between the two modes is weak, i.e., 7 < 7c = 
min(2Ka//3, 2r][3/a), one has a solution that both ur and 
Uip are non-zero as schmeatically shown in Fig. |3] (a). 
However, when 7 > 7c, one has ur=0 and u^ = — /3/2k, 
if /3/2k < a/2ri, (u<^=0 and u^ = -a/2r], if I3/2k > 
a/277), i.e., one of the soft mode is fully suppressed by 
the other mode that has stronger instability, because of 
the coupling between the two modes, as shown in Fig. [3] 
(b). This is exactly the case in CaMnOa. In bulk mate- 
rials, the AFD mode is more unstable, therefore the FE 
mode is fully suppressed. In the 1:1 SL, when the Mn02 
intra-layer spin configuration is AFM, the AFD insta- 
bility a is still larger than that of the FE mode despite 
of the interface effects, therefore, there are still signifi- 
cant AFD rotations, even though they are smaller than 
those in bulk CMO. However, when the Mn02 intra-layer 
spin confi gura tion changes to FM, the FE instability is 
enhanced [20[, and becomes stronger than the AFD in- 
stability with the help of the BaO layer pinning effects 
to the AFD mode. Therefore, the AFD mode is fully 
suppressed. The suppression of the AFD mode greatly 
enhances the FE modes, again because of the coupling 
between the FE and AFD modes. The results for the 2:2 
SL can also be understood in the same scenario. 

The magnetoelectric coupling effects in the 
CMO/BTO SLs can be observed by various experi- 
mental techniques. For example, one should be able to 
observe large polarization change by applying magnetic 
field. This can be best seen near the paramagnetic (FM) 
to AFM phase transition temperatures, where relative 
small magnetic field is needed. Or one may simply ob- 
serve the polarization change at FM to AFM transitions. 



Beside the electric polarization, the dielectric constants 
are also expected to change dramatically under magnetic 
field near the magnetic phase transitions. 

To summarize, we have demonstrated a novel mecha- 
nism that could lead to giant magnetoferroelectric cou- 
pling in the multiferroic CMnOs/BaTiOa superlattices. 
The key idea is that the energetics of the instabilities, 
such as antiferro-distortive mode, ferroelectric mode, and 
the magnetic ordering can be drastically modified by in- 
terface engineering to enhance the magnetoelectric cou- 
pling. The enhancement of the magnetoelectric coupling 
is due to the interface atomistic effects which could be 
much stronger than those with mechanical coupling. It 
therefore opens a new path to design novel multiferroic 
superlattices with strong magnetoelectric coupling. 
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